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ABSTRACT 

Recent results on the non-universality of the Initial Mass Function (IMF) have shown 
strong evidence of IMF variations with galaxy velocity dispersion, with a correspond- 
ing impact on other stellar population parameters, line indices and colours. Using a set 
of stellar population models with varying IMF slopes, we provide additional caveats 
on the assumption of a universal IMF. The present study shows that the derived star 
formation histories of early- type galaxies vary significantly with the IMF slope. For 
instance, a steepening in the slope of a single power-law IMF decreases substantially, 
by a factor of up to four, the contribution of the old stellar populations to the total 
light/mass. This trend is milder for a segmented-like IMF shape, where the contribu- 
tion of the very low mass stars is decreased. It is also shown that, by tuning each IMF 
slope to its prescribed value according to each galaxy velocity dispersion, a sample of 
early-type galaxies covering a range of masses yield comparable star formation histo- 
ries. On the one hand, a small contribution from relatively young stellar populations 
appear in the star formation histories of most massive elliptical galaxies when adopt- 
ing a steep IMF. On the other hand, we find that low mass early-type galaxies that 
look like genuinely young objects with a standard IMF (i.e. ^^baby elliptical galaxies ^^), 
turn out to be older when a slightly flatter IMF is employed. In summary, the use of 
a non-universal IMF, tuned according to the velocity dispersion of the galaxy, seems 
to provide more consistent results. 
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1 INTRODUCTION 

The Initial Mass Function (IMF), which is the distribution 
of stellar masses in a single population at the time of 
birth, has been usually assumed to be universal. This 
means that we consider the same IMF at all cosmic times, 
no matter the environment or the mass of the galaxy. 
This universal assumption was adopted mainly due to the 
difficulties in obtaining direct constrains. So far, the IMF 
can only be directly constrain ed from our ow n Milky Way 
by counting individual stars (Sal peterl ll955l V Many works 
have supported this universality, warning us about the 
consequences of freely varying i t (see e.g. iGilmord 
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Narayanan &; Pavel 12013 ). For example, an IMF with a 



steeper slope than Salpeter would probably fit better the 
near-infrared Call triplet, but it would also produc e too red 
V-K colours and too high M/L fe.g. ICenarro et al.ll2003h . 
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Elbaz, Arnaud & Vang;ioni-Flam (1995), Vazde kis et al.l 
()l996i ) or IVazdekis et al.l (|l997l ). among others. Re- 
cently, abundant observational evidence from a vari- 
ety of independent techniques have pointed out that 
a non-universal IMF might be in fact the real case 
(|Tortora, Romanowskv Napolitanol |2012| and references 
therein). From the measurements of spectral features sensi- 
tive to the slope of the IMF in massive giant elliptical galax- 
ies (e.g. Call triplet, the Wing- Ford band and the NaI8200), 
it has been discovered the existence of a larger population 
of low-mass stars th an predicted by a Kroupa or even 
a Sa l peter IMF (e.g. IVazdekis et al.l l2003l. ICenarro et al " 



20031 . iFalcon-Barroso et al.l l2003l. Ivan Dokkum &: Conrov 
2011 IConrov van Dokkuml l2012al (CvD12 he re after) 
Spiniello et aD l20ll Ivan Dokkum Conrovl l2012l . 
Conrov van DokkumI l2012bl ). Colours are also highly 
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sensitive to IMF variations, in the sense that redder galax- 
ies have bottom-heavier IMFs, i .e. steeper slopes (see e.^. 
Dutton. Mendel Simard' 2012, Pforr. Maraston Toninil 
2012 . Ricciardelh et al. 2012, Vazdekis et al. 20li) In 
addition, recent works have added strong evidence on 
the non- universality of the IMF by constraining the dark 
matter fraction in galaxies, with gravitational lensing 
and dynamical studies, f inding an IMF dependenc e with 
galaxy mass fe.g: Grillo et alJ l2009l. lAuger et al.l | 2010| . 



Napolitano. Romanowskv TortoraEoiOljTreu et al.l 201' 
Barnabe et al. l201lL Thomas et alT201lL ICappellari et al l 
2012bl. button et all l20ll ISonnenfeld et al.l l2012l 



Tortora, Romanowskv &: N apolitanol |2012| V Moreover, in 
Ferreras et al 1 (|2013l ) (F13 hereafter) the authors have 



shown a tight correlation between the velocity dispersion 
of early- type galaxies and the IMF slope, in the sense that 
higher velocity dispersion galaxies need ste eper IMF slopes, 
while t he low mass ones need flatter slopes. ICappellari et al.1 
(|2012al ) have shown a similar result, where the use of dif- 
ferent IMF slopes varies the amount of dark matter and 
implies a change in the dynamical masses by a factor of 
three in early-type galaxies. 

As the IMF is a parameter that clearly influences most of 
the stellar population properties, its determination is crucial 
to understand the formation and evolution of galaxies. Most 
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many others). However, a number of still open questions 
could be addressed by changing the IMF according to the 
recent claims in this direction. 

The work here conducted has been devised as a " user-guide " 
for observers, who usually apply their analysis tools using 
models a s sumin g a standard, universal I MF s uch as Kroupa 
( Kroupal l ioOlh , Chabrier (jChabried l2003l ) or Salpeter 
( SalpeteJll955h . Rather than prescribing a particular IMF, 
it is shown here that more consistent results can be obtained 
when the IMF slopes are tuned according to the velocity 
dispersion of the galaxies. Following the trend in F13, we 
will use throughout this work the notation /i(cr) for these 
cr-dependent IMF slopes. We analyze the behavior of the 
Star Formation History (SFH) and the derived stellar mass 
(Mstar) of Early- Type Galaxies (ETGs) when models with 
different IMF slopes and shapes are used. 
Section 2 describes the high quality long-slit data sample 
used to derive accurate SFHs for galaxies covering a wide 
range of masses. In Section 3 we analyze the impact of 
varying the IMF slope on various galaxy properties such 
as the M/L ratio, the SFH, the derived ages or the stellar 
masses. Section 4 is dedicated to an illustrative case and our 
conclusions and summary are found in Section 5. We have 
adopted a concordance cosmology with r2rn= 0.3, 11a=0.7 
and Ho=70kms"^ Mpc~^ 



2 THE DATA 

Spectra from two ETGs samples are used for the present 
study. The main sample, with high-quality long-slit spectra 



make the bulk of the work, while a second sample with fiber 
spectra from the SDSS {Sloan Digital Sky Survey) is used 
for statistical purposes. 

The first sample , with S/N^40A~^, was c onstructed 
with spectra from | Sanchez-Blazquez et al. (l2006l) (P SB06) , 
lYamada et al. 1 (|2006l ) (YOG) and iFerre-Mateu et al.l (|2012h 
(AFM12). The reader is referred to these papers for a com- 
plete description of the data, their quality and reduction 
process. The first two references contain ETGs covering a 
range of masses (cr= [50-300] kms""*^), while the third one in- 
cludes massive compact galaxies that are found in the Local 
Uni verse. The la t ter ar e a family of unique galaxies found 
in iTruiillo et al.1 (120091) and characterized in AFM12 and 
iTruiillo, Carrasco &: Ferre- Mateul (|2012l ). All these galaxies 
allow us to cover an ample range of properties, like masses, 
luminosities, sizes, stellar populations, etc. Therefore, this 
sample is ideal to detect the effects of IMF variations. In 
particular, the YOG sample covers the colour-magnitude di- 
agram of the Virgo Cluster. At the very high mass regime 
(cr ^ 300 kms~^) galaxies with the steepest IMF are located, 
according to F13 and CvD12. The local compact massive 
galaxies were selected because they have revealed unusually 
large fractions of young components on their derived SFHs. 
The massive ellipticals were selected to have similar veloc- 
ity dispersions to the compact galaxies, albeit larger radii, 
to compare their mass estimates (cr= 200 km s~^). The low- 
mass ellipticals are useful probes to study the SFH differ- 
ences, as a fraction of them show young ages, comparable to 
those in the massive compacts. Table 1 lists all the galaxies 
in the high-quality spectra sample with their main proper- 
ties. Despite their different origins, they were all treated and 
analyzed in a similar way. 

The sample of ETGs from t he SPIDER projec t with SDSS 
spectra of 40,000 objects (La Barbera e t al.l l2010l , LBIO 
hereafter) , has been used as the benchmark for sanity checks 
performed to test the results obtained. It was also used to 
find rather extreme, but illustrative, cases in Section 4. 



3 GALAXY PROPERTIES WITH VARYING 
IMF 

A study on the quantitative impact of varying the IMF shape 
and slope on the stellar population properties is presented in 
this section. A first consideration should be done regarding 
the shape of the IMF assumed for the models. Historically, 
two alternative IMFs have been considered: a single-power 
law ove r all stellar ma sses, characterized by its log arithmic 
slope /i (Salpeter"l955') and a multi-segmented IMF (Kroupa| 
2001). These two IMFs will be regarded as standard through 
this work. For our illustrative purpos es, th e other popular 
choice, the log-normal IMF ( Chabrierll2003h is very similar 
to that of Kroupa and will not be considered here. 



3.1 Mass-to-light ratios 

We can study the dependence of the mass-to-light ratio 
(M/L) with the IMF using stellar population synthesis mod- 
els. We have selected three representative ages (1, 2 and 
lOGyr) and three metallicities ([Z/H]= -0.40, 0.00 and 
0.22) from the extended M ILES single-burst models (MIUS- 
CAT: IIVazdekis et al.llioTj , V12 hereafter). We have selected 
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Table 1. Galaxy sample of high-quality long-slit spectra 



Galaxy ID Source environment type z cr (km s ^) 



NYU 54829 


AFM12 


Field 


MC 


0.085 


137 


NYU 321479 


AFM12 


Field 


MC 


0.128 


221 


NYU 685469 


AFM12 


Field 


MC 


0.149 


204 


NYU 796740 


AFM12 


Field 


MC 


0.182 


203 


NYU 890167 


AFM12 


Field 


MC 


0.143 


233 


NYU 896687 


AFM12 


Field 


MC 


0.130 


223 


NYU 2434587 


AFM12 


Field 


MC 


0.172 


206 


NGC4365 


Y06 


Virgo 


HME 


0.0041 


265 


NGC 4472 


Y06 


Virgo 


HME 


0.0033 


300 


NGC 2329 


PSB06 


Coma 


ME 


0.0192 


225 


NGC 4473 


Y06 


Virgo 


ME 


0.0074 


193 


NGC 4621 


PSB06 


Virgo 


ME 


0.0013 


230 


NGC 4697 


Y06 


Virgo 


ME 


0.0041 


180 


NGC 5812 


PSB06 


Field 


ME 


0.0065 


215 


NGC 4239 


Y06 


Virgo 


LME 


0.0031 


63 


NGC 4339 


Y06 


Virgo 


LME 


0.0043 


114 


NGC 4387 


Y06 


Virgo 


LME 


0.0015 


98 


NGC 4458 


Y06 


Virgo 


LME 


0.0210 


105 


NGC 4464 


Y06 


Virgo 


LME 


0.0041 


135 


NGC 4467 


Y06 


Virgo 


LME 


0.0047 


69 


NGC 4489 


Y06 


Virgo 


LME 


0.0032 


52 


NGC 4551 


Y06 


Virgo 


LME 


0.0039 


105 



Main galaxy propertie s: (1) galaxy I D; (2) publication source: PSB06 (ISanchez-Blazauez et al.ll2006l ). Y06 (jYamada et al.ll2006l ) and 

AFM12 ( Ferre-M ateu et al.|[20l3 ): (3) environment of the galaxy; (4) galaxy type: MC (Massive Compact), HME (High Mass 
Elliptical), ME (Massive Elliptical) or LME (Low-Mass Elliptical); (5) redshift; and (6) central velocity dispersion. Values are taken 
from the publication source or from the LEDA database ( http: / /leda.univ-lyonl.frT]) when not available. 



six different IMF slopes: /i = [0.3, 0.8, 1.0, 1.3, 1.8 and 2.3] 
both for a single- and segmented-p ower law (unim o dal an d 
bimodal following the notation of Vazdeki s et al.j (|l996')). 
As a reference, /i= 1.3 corresponds to a Salpeter IMF slope 
with a unimodal shape. Instead, a bimodal IMF with the 
same slope resembles closely to a Kroupa IMF. In addition, 
a unimodal IMF with slope /x ^1.0 would render similar 
results in the studied properties to the Chabrier IMF. 
Figure 1 shows the dependence of the M/L with the slope of 
the IMF for stellar populations of different ages with both 
shapes. It can be seen that the M/L from the stellar pop- 
ulation models presents a minimum. For the M/L calcula- 
tion both the alive stars contributing to the light and the 
stellar remnants are taken into account. Whereas the steep- 
est IMFs emphasize the contribution of long-lived low-mass 
stars, the flatter IMFs favour the number of remnants re- 
sulting from short-lived massive stars. Note that if the stel- 
lar remnants were not included, the resulting M/L would 
not show a minimum when varying the IMF slope (see the 
predictions shown in the MILES webpag43)- The location 
of this minimum depends on the age of the SSP: younger 
ages reach their minimum at flatter slopes. In the unimodal 
case, the minimum is found at /x ^ 0.8 for young popula- 
tions, whereas for it is located at fi^l.O for the old ones. In 
the bimodal case the minimum is seen at /x^l.O (1.3) for 
the young (old) populations. It is worth noticing the strong 
variation of the M/L with the IMF slope. For the unimodal 
shape, the M/L value increases strongly for both the flattest 
and the steepest IMFs. For the bimodal IMF a significantly 
milder variation is seen, particularly for the oldest stellar 

^ http://miles.iac.es/ 



populations (see also F13). The choice of one shape over the 
other will be affecting in a different way the derived stel- 
lar populations properties, such as the M/L. The plot also 
shows a trend with the broad band filters. At a fixed IMF 
slope, lower M/L values are obtained for the redder filters. 

3.2 Star Formation Histories 

To obtain the SFHs of our high-quality galaxy sample we 
apply a ful Lspectrum- fitting approac h. We use the code 
STARLIGHT (|Cid Fernandes et al.ll2005l ) with the MIUSCAT 
SSP SEDs. This code fits a set of SSP models to the 
observed spectra and selects the best population mixture. 
The SFH is derived for each galaxy using different base 
models with all the IMF slopes selected above, both for the 
unimodal and the bimodal cases. 

The high quality of our spectra allows us to perform an ac- 
curate fitting. We are aware that our results may be biased 
due to the relatively short spectral range of the spectra 
(3800-5300 A). Although this range encompasses most of 
the commonly used absorption line features, the reported 
IMF-sensitive indices are located at redder wavelengths. 
However, the blue part of the spectrum also shows some 
sensitivity to the IMF, as shown in figure 9 of VI 2. To 
check the robustness of our results, we have performed 
various tests to take into account different sources that 
might change our results (see Appendix A). We first masked 
different features in the spectra, finding that the variations 
on the mean derived ages are negligible (smaller than 2%). 
We also compared the results from selecting a short versus 
a longer spectral range, which may account for variations 
in the derived luminosity(mass)-weighted ages of the order 
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Figure 1. The M/L ratio in several broad-band filters is plotted against the IMF slope for models of 1 Gyr (black), 2 Gyr (blue) and 
10 Gyr (red) and for different metallicities ([Z/H]= -0.40, dotted fine; [Z/H]=0.00, dashed-dotted fine and [Z/H]=0.22, dashed fine). 
Upper rows correspond to models with a unimodal IMF and lower ones to models with a bimodal IMF. A minimum on the M/L is seen 
around ^ ~ 0.8-1.0 (unimodal) or ^ ~ 1.0- 1.3 (bimodal), depending on the SSP age. 



of a 8% (12%). Finally, due to the se nsitivity of some 
colou rs to the IMF variations (see e.g. iRicciardelli et al.l 
I2OI2I I. we have compared the synthetic colours obtained 
from our derived SFHs to those from the SDSS broadband 
photometry. Note that the latter extend to the UV and the 
IR, i.e. sampling the different populations in the galaxy. We 
find a good agreement between both. These tests confirm 
the robustness of the results and trends obtained in Section 
3. From now on, all the SFHs presented here are derived 
using the spectral range of our high quality long-slit sample. 

For illustration. Figure 2 shows the SFHs of four galaxies, 
colour-coded by galaxy family: the very massive galaxy 
NGC4472 (black), the massive elliptical NGC4473 (green), 
the massive compact NYU 2434587 (purple), and the 
low-mass galaxy NGC4458 (yellow). The panels show 
the derived SFHs with base models where the IMF slope 
steepens from left to right. In the upper rows we show 
the results for the unimodal case and in the lower ones 
for the bimodal one. The dashed vertical line shows the 
mean mass- weighted age for each varying slope, while the 
solid line is the mean value derived from the standard 
/i= 1.3 slope. The SFHs for the other galaxies are shown in 
Appendix B. 

Some interesting conclusions are drawn from Fig. 2. First, 
we see noticeable differences in the derived SFHs for these 
types of galaxies: while massive ellipticals are mainly old, 
both in light and mass, massive compacts show an unusually 
large fraction of young stellar populations (^2 Gyr), both 



light- and mass- weighted. We refer the reader to AFM12 
for an extensive discussion on this issue. Second, we see 
a general trend of decreasing the fraction of old stellar 
populations with increasing IMF slope. This trend is shown 
in Fig. 3 (upper panel), where we indicate how the fraction 
of stellar populations larger than 5 Gyr evolves with the 
IMF slope. 

This trend is seen in all galaxies and for both IMF shapes 
studied here, though it is emphasized for the unimodal 
case. For various massive galaxies a stellar population with 
age ^5 Gyr emerges when adopting an IMF steeper than 
Salpeter, leading to a more complex SFH. The impact is 
significantly milder for the bimodal case. Note however, that 
according to F13 the best fits to the various IMF-sensitive 
features for the most massive galaxies require an IMF 
slope around 2 and 3 for the unimodal and bimodal shape, 
respectively. 

This rejuvenating effect implies a decrease in the mean 
age of the galaxy, as shown in the lower panel of Fig. 3, 
where the derived mean mass- weighted ages are shown with 
respect to the IMF slope. It is seen that in the unimodal 
case, the derived ages might decrease by a factor of as much 
as 3 from ^— 0.8 to \i— 2.3, while for the bimodal case is 1.5. 

As stated before, recent results have claimed that more 
mass ive galaxies demand steeper IMF slopes than Salpeter 
(e.g. ICappellari et al.ll2012al. Ivan Dokkum h Conrovll20ll 
CvD12. ISpiniello et al.1 12012 ). including a systematic cor- 
relation with the found by F13. These papers pose 
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an upper limit for the steepest slope around /i ^2 for a 
single power-law IMF shape. Therefore, for our illustrative 
purpose we will only discuss the unimodal case in this 
part of the work. Which SFHs will be derived for our 
objects if we follow the /x-cr relation? We find that most 
massive galaxies for which we adopt /i(cr) ^2 are very 
old. However, for most cases, a non-negligible fraction of 
relatively younger stellar components tends to appear (up 
to 10%). This also applies to the massive galaxies, for 
which a Salpeter or slightly steeper IMF is demanded. On 
the contrary, the contribution from old stellar populations 
increases for the low-mass ETGs if adopting flatter slopes 
than Salpeter (/i -0.8-1.0). 

In summary, our results suggest that the SFHs of ETGs 
of all masses tend to converge to a common pattern, when 
the prescribed IMF slopes, according to the galaxy velocity 
dispersion, are used. This common SFH pattern involves 
a varying amount of recent residual star formation. The 
SFHs of ETGs have always been a challenge in modern 
astronomy, as they represent a direct way to constrain 
competing galaxy formation models. Although it is well 
established that the bulk of the stars in ETGs formed at 
high redshift (z'-^2), a small "frosting" of youn^ stars forni s 
on top of it at more recent epochs (e.g: iTrager et al. I l2000h . 
Observations in the UV are particularly useful to trace these 
recent bursts, as the rest-fra me UV is mostly sensitive to 
stars l ess than —1 Gyr old (e.^:lEllis, Abraham Dickinson 



2001 . iKaban. Nolan Ravchaudhurvl l2005l . iKavirai et aL. 
20051 . lYi et al.ll2005l . iFerreras et al.ll2006l ). All these works 
show that there is inde ed this residual star fo r mation, which 
evolv es with redshift (jFukugita et al.l |2004| . iKavirai et al.l 
20081). It can be up to —10% in the local Universe (e. g 
Kavirai et alll2007l . iNolan. Ravchaudhurv KabanI [ioOTh . 
in agreement with our results. 

Nonetheless, for the compact massive galaxies with an 
unusually large contribution of young stellar populations, 
as reported in AFM12, the derived SFHs do not change 
substantially if the IMF is varied according to the /i-cr rela- 
tion. This also applies when adopting the more conservative 
bimodal IMF shape. Therefore, these galaxies further show 
their unique nature in the present analysis. 
Table 2 summarizes all the mass-weighted properties 
derived from the SFHs that are obtained when adopting 
for each galaxy their /i(cr). The other derived quantities 
for both the two IMF shapes adopted in this work can be 
found in Table Bl and B2 in Appendix B. 

3.3 Stellar Masses 

Stellar masses are computed from the mass-to-light rela- 
tion corresponding to the SSP mixtures that match best the 
galaxy spectra, as obtained from the full-spectrum-fitting 
technique applied here. For our massive compact galax- 
ies, the M/L values were calculated in t he r-band of the 
SPSS. We adopt ed the GALFIT magnitude (|Peng et al.ll200l 
IPeng et al]|2010l l as luminosity estimate from AFM12, after 
correcting for galactic extinction. For the elliptical galax- 
ies we adopted the M/L values in the V band as lumi- 
nosity estimate. We used the values tabulated for the to- 
tal magnitude, corrected in the same band, from LEDA 
(http:/ /leda.univ-lyonl.fr/). Figure 5 shows the dependence 
of the derived stellar mass with the IMF slope. This is the 



observable with the largest sensitivity to the assumed IMF 
shape and slope. It is seen that for the unimodal IMF the 
stellar mass grows up by a factor of 3 from the standard 
value /i= 1.3 to the steepest case /i=2.3, with the smallest 
stellar mass obtained at /i ^ 1.0. For the bimodal IMF shape 
the impact is significantly milder. It shows a smaller increase 
of the factor of 1.5, emphasized for the massive galaxies. The 
trend of increasing mass with a steepening of t he IMF slope 
is analog to the results found in Cappell ari et al.l (l2012al) 
for the dynamical masses. This finding provides further sup- 
port for the need of using the appropriate IMF, in order to 
avoid galaxy misclassifications due to wrong stellar masses. 
As an example, assume that a sample of massive galaxies will 
be constructed under the condition: M* > 10^^ ^0- If 
stellar population parameters were derived with a Chabrier 
IMF, the stellar masses would be underestimated and a frac- 
tion of galaxies would be lost. 



4 NGC4489: BABY GALAXY OR 
NON-STANDARD IMF? 

We have seen in Section 3.2 that by tuning the IMF slope 
related to the velocity dispersion of the galaxy, according 
to recent evidence in this direction, we tend to derive more 
similar SFHs among galaxies of different masses. However 
such a choice does not change drastically the results 
one would infer from adopting a standard-universal IMF 
(e.g. Salpeter or Kroupa). We find that the most massive 
galaxies are nevertheless old, despite IMF- modulated young 
population contributions. However, here we will show that 
there are galaxies for which the use of a standard IMF 
would lead to surprising results. 

Here we focus on the low-mass elliptical NGC4489. This 
galaxy has been previously reported to have young mean 
luminosity- weighted ages from l ine-strength studies (^3- 

5 Gyr; e.g: iGorgas et al.1 Il997l. iTerlevich Forbed pQOi . 



ICaldwell. Rose Concannonll2003l . iTantalo ChiosilboOJ . 
YOG). However, this does not necessarily imply that the 
galaxy is composed of genuinely young stellar populations, 
as the old component can be masked by the luminosity of a 
tiny fraction of young stars on top of it. The mass- weighted 
SFH reveals the composition of NGC4489, as presented in 
Figure 5. It shows, like in Fig. 2, the derived SFHs with 
different IMF slopes. From the Salpeter panel (unimodal 
/i= 1.3) one would infer that NGC4489 is a genuinely 
young object with almost no contribution from stellar 
populations older than 5 Gyr (less than 5%), i.e. a ''baby 
elliptical galaxy". However, considering a slightly flatter 
IMF slope as claimed fo r an elliptical galaxy of such small 
velocity dispersion fe.g. ICappellari et al.l bOlial . F13), an 
old component appears, which contributes ^ 65% to the 
total mass. This result is in much better agreement with 
those obtained in Section 3.2. 

Do other galaxies like NGC 4489 exist in the Local Universe 
that can be considered as baby elliptical galaxies assuming a 
standard IMF slope? Out of the SDSS-based LB 10 sample, 
we selected those galaxies with a mean mass-weighted 
age younger than ^ 5 Gyr (Salpeter IMF). Only 11 baby 
elliptical galaxies were found, with a similar behavior to 
NGC 4489, as shown in Fig. 6. Once again, the apparent 
extreme youth of these galaxies results from the inadequate 
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Figure 2. Star Formation Histories for the very massive galaxy NGC 4472, the massive eUiptical NGC 4473, the massive compact galaxy 
NYU 685469 and for the low-mass elliptical NGC 4458. From left to right, the mass-weighted SFHs derived with increasing IMF slope 
(/i), as indicated on the top of each panel, for the unimodal in the upper row and the bimodal in the lower one. The dashed vertical 
line corresponds to the mean mass-weighted age derived for the galaxy with each slope, while the solid line represents the mean value 
derived from the standard ii= 1.3 slope. The fraction of mass-weighted old populations (>5 Gyr) is quoted in each panel. Note the trend 
of a decreasing old stellar population contribution when steepening the IMF slope, more remarkably for the massive compacts and the 
low- mass galaxies. 
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Figure 5. Derived SFHs for NGC 4489 for increasing IMF slopes, 
as quoted within each panel as in Figure 2. Note that for a stan- 
dard Salpeter IMF (third panel) this object is genuinely young, 
while a slight flattening of the IMF slope gives a SFH with a large 
contribution of old stellar populations. 



use of a standard-universal IMF. 



5 CONCLUSIONS 

We have gone one steep further on the discussion about 
the universality of the IMF. Although in this work we do 
not attempt to constrain the IMF, we explore and quantify 
the impact of varying its slope and shape on various galaxy 
properties (M/L, SFHs, Mstar and mean ages). These are 
summarized as follows: 

• The choice of the IMF shape (unimodal vs bimodal) 
results in different stellar population parameter estimates, 
in particular for the steepest slopes, as seen from the M/L 
ratios. The unimodal (bimodal) IMF shape shows stronger 
(milder) deviations on the derived properties with respect 
to the standard IMF. 

• We find a minimum on the M/L for an IMF slope that 
depends slightly on the age of the SSP: /i = 0.8 for young 
SSPs, increasing towards /i= 1.0 for the oldest stellar popu- 
lations in the unimodal case. For the bimodal this minimum 
is located at /i=1.0. This implies that at /x 0.8-1.0, the 
stellar mass derived from the SFHs will be minimized. In 
the unimodal case the use of steeper IMF slopes increases 
the derived stellar masses by a factor of 4 (from /i= 0.8 to 
/i=2.3), while the increase is by a factor of 1.5 in the bi- 
modal one. This is the most sensitive parameter to the IMF 
variations. 
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Figure 3. Fraction of old populations (>5 Gyr) and mean ages derived from the mass- weighted SFHs using different IMF slopes. Symbols 
follow the previous colour code: high massive galaxies in black, massive ellipticals in green, the massive compacts in purple and the low- 
mass ellipticals in yellow. The vertical solid gray line indicates the position of the standard slope for both IMF shapes {12= 1.3). The 
bimodal case shows milder variations with a decrease of a factor of 1.5, half the variation observed for the unimodal case. 
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Figure 4. Stellar masses obtained from the derived M/L ratios with different IMF slopes. Note the large difference between unimodal 
and bimodal shapes, in particular for the steepest slopes. The stellar masses increase by a factor of 4 in the unimodal case, while it only 
accounts by a factor of 1.5 in the bimodal case. 



• IMF slopes steeper than Salpeter tend to produce SFHs 
with smaller contributions of old populations. This effect 
decreases the mean age of the galaxy, both in mass and 
light. This effect is observed for both IMF shapes analyzed 
here, although in the unimodal case the decreases are more 
drastic, by a factor of up to 3 (from /i=0.8 to /i=2.3). It 
shows milder versions, up to 1.5, for the bimodal case. 

• The choice of an IMF slope according to the velocity 
dispersion of the galaxy seems to provide more constrained 
results. While standard IMFs produce rather mass (and age) 
dependent SFH types, the use of the /i-cr relation produces 
comparable SFHs for the ETG family. The most massive el- 
lipticals, for which we adopt slopes steeper than Salpeter, 
present a small, but not negligible contribution of relatively 
younger stellar populations. Within this scenario, massive 
galaxies would no longer be completely dead and passively 



evolving old objects, supporting recent claims in this direc- 
tion. Adopting flatter IMF slopes for the low-mass galaxies 
will render more similar, but still different, SFHs to the mas- 
sive ones. However, the SFHs derived for our compact mas- 
sive galaxies, which show unprecedented large contributions 
of young stellar populations, do not change significantly if 
the IMF slope is varied. We also find galaxies with appar- 
ently no contribution from old stellar populations using a 
Salpeter IMF {hahy elliptical galaxies). However, the use of 
flatter slopes, implied by their low velocity dispersion, un- 
cover a conspicuous old population. 

These results point out that a variation of the IMF 
slope might solve, or at least alleviate, some of the problems 
encountered when analyzing galaxies of varied masses. 
Contrary to earlier warnings against the use of non stan- 
dard IMFs, the present study provides new caveats on the 
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Table 2. Stellar properties with /j.(a) 
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Stellar populations parameters estimated from the derived SFHs using the prescribed IMF slope, according the the galaxy velocity 
dispersion (/i(cr)): (1) galaxy ID; (2-3) The selected /i((T) for each galaxy, both for the unimodal and the bimodal IMF. We consider the 

closest slope in our model (0.8, 1.0, 1.3, 1.8 or 2.3) compared to the value obtained following F13 equation (quoted within the 
parenthesis). Note that all the low-mass galaxies have /i(cr)= 0.8, because F13 equation is valid for galaxies with a >150 kms"-*^; (4-5) 
Mean mass- weighted ages; (6-7) Fraction of old stellar populations (>5Gyr); (8-9) Stellar mass derived from the SSP mixture resulting 

from the derived SFH. 



assumption of a universal IMF, such as the questionable ex- 
istence of early-type galaxies with virtually no contributions 
of old stellar populations, or the misclassified non-massive 
galaxies due to mass underestimation. Moreover, the 
complex SFHs found for most massive galaxies when using 
steeper IMFs open new means to explain their formation 
and evolution. 
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APPENDIX A: TESTING OUR METHOD FOR 
DERIVING THE SFHS WITH STARLIGHT 

Full-spectrum fitting codes are known to be slightly sensitive 
to the initial conditions, such as the spectral range (e.g. blue 
vs red), or masked features (e.g. potential emission lines). 
We have tested the robustness of our results by varying these 
initial conditions, while keeping the IMF shape and slope 
fixed to a unimodal /i=1.3. 

Altogether, these tests reveal that the trends and estimates 
derived from the relatively short spectral range of the long- 
slit galaxy spectra are robust. 



Al Masking test 

Three masking configurations were used, involving different 
features on the spectra of four of our galaxies, as illustrated 
in the upper panel of figure Al: (a) orange - no mask is used; 
(b) blue - we mask all the features in the spectra, both po- 
tential emission and Balmer-absorption lines; (c) green - we 
only mask the Balmer-absorption lines. The resulting mean 
luminosity-weighted ages could be used to provide a more 
quantitative estimate of the uncertainties derived from the 
method. We find that these mean ages, which are indicated 
for each panel, vary insignificantly, by less than 2%. 
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Figure A3, g-r and r-i colours derived from our SFHs on the 
short spectral range (open green circles) compared to those re- 
ported by the SDSS (blue filled circles). 

this filter (see iRicciardelh et al.l 120121 ). 



A2 Short vs Long spectral range 

We have compared the mean ages derived from the relatively 
short range (3800-5300 A) of the high quality long-slit spec- 
tra of our sample of galaxies with those derived from the 
full SDSS spectral range. The latter covers bluer and redder 
wavelengths (3500 to 8500 A ), better sampling the differ- 
ent contributions from the various stellar populations. For 
this test we have selected a set of ETGs from the sample of 
LB 10, covering a wide range in ages. Moreover, as the com- 
pact massive ga laxies were initiall y selected from the DR6 
SDSS database (jTruiillo et al.|[2009l l. we also compared their 
results. The luminosity- weighted estimates show a rather 
good agreement, with deviations on the derived mean age 
smaller than 8%, whereas the mass- weighted estimates show 
deviations on the mean age of around 12%. The larger vari- 
ations are found for the galaxies with the oldest ages, as 
expected for the errors associated to this age regime. 



APPENDIX B: SFHS DERIVED FOR THE 
SAMPLE OF HIGH-QUALITY SPECTRA 

In this Section we show the SFHs derived for all the galax- 
ies of our sample with high-quality spectra. We do not in- 
clude here the galaxies already shown in Fig. 2 and 5. The 
first six figures correspond to the SFHs of the massive com- 
pact galaxies. The next six show the SFHs of the massive 
and very massive ellipticals and the last six correspond to 
the low-mass ellipticals. As in Fig. 2, the panels show, from 
left to right, the derived SFHs obtained when using models 
with increasing IMF slope. In the top panels we show the 
results when adopting a unimodal IMF, whereas in the bot- 
tom panels we show the bimodal case. The dashed vertical 
line corresponds to the mean mass- weighted age derived for 
the galaxy with each slope, while the solid line represents 
the mean value derived from the standard ii— 1.3 slope. 



A3 Broadband photometry vs Spectrum fitting 

Colour variations with the IMF slope ha ve been 

already reported (e .g. Ricciardelli et al.l |2012| , 

IPforr. Maraston k Toninil I20121 K However, the effect 
of the colours in our short spectral range is almost negligi- 
ble. We compared the SDSS photometric colours with those 
obtained from the SFHs derived from the short spectral 
range of these galaxies in LB 10. To obtain the synthetic 
colours we use the user-friendly web-tool facility dedicated 
for this purpose in the MILES web-page "Get spectra from 
a SFH" . This facility provides us with the indices, masses 
and colours for our derived SFHs. Figure A3 shows that the 
agreement is remarkable between the derived and observed 
g-r and r-i colours. The web-tool does not provide us with 
the magnitude in the u band, because the MIUSCAT SEDs 
do not cover the whole spectral range required to measure 
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Figure Al. Different masking configurations employed for the spectrum- fitting are indicated on the upper panel. The vertical blue and 
green dashed lines mark the masked lines for each configuration. The SFHs derived for four elliptical galaxies are shown in the following 
rows of panels. No masking (left), masking emission and Balmer lines (center) and only masking the Balmer lines (right). The mean 
mass-weighted age is shown for each SFH, giving us a quantitative estimate of the uncertainties of this method. 
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Figure A2. Mean mass- and light- weighted ages (left and right, respectively) obtained from the SFHs derived using the short and large 
spectral ranges. Filled symbols correspond to the galaxies from the LB 10 sample, while open symbols correspond to the compact massive 
galaxies with SDSS data. The vertical lines on the spectra indicate these two spectral ranges. 
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Figure Bl. Star formation histories derived for the galaxies in our high-quahty sample. The IMF slope of the adopted base models 
increase from left to right as quoted on each plot, for the unimodal in the upper row and for the bimodal in the lower one. The dashed 
vertical line indicates the derived mean mass-weighted age with each slope, while the solid line represents the mean value obtained from 
the standard IMF slope (1= 1.3 

. As for the previous plots, the SFH is colour-coded by galaxy type: black for the high mass ellipticals, green for the massive ellipticals, 
purple for the compact massive galaxies and yellow for the low-mass ellipticals. A trend where the contribution of the old stellar 
populations decrease with increasing IMF slope is clearly seen for all the galaxies. 
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Figure B2. Fig. Bl continued 
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Figure B3. Fig. Bl continued 
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Table Bl. Variations on the stellar population parameters when adopting a UNIMODAL IMF 



Age (Gyr) %old Mass (x IQ^^Mq) 

Galaxy ID IMF Slope IMF slope IMF slope 
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1.0 


1.3 


1.8 


2.3 
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1.8 


2.3 
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1.0 


1.3 
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2.3 
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63 
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14.2 


13.7 
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7.8 


90 
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5.6 


3.4 


60 


55 


50 


30 


20 


1.03 
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55 
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9.1 
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91 


86 
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9.9 


9.7 
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5.0 
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78 


60 


35 
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3.39 


5.79 


NGC 5812 
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3.1 
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91 
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71 
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1.46 


NGC 4387 


13.7 


12.9 
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Mass- weighted age, fraction of old stellar populations (>5 Gyr) and stellar mass obtained from the derived SFH with a unimodal IMF 

shape. 
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Table B2. Variations on the stellar population parameters when adopting a BIMODAL IMF 
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Age (Gyr) %old Mass (x IO^IMq) 

Galaxy ID IMF Slope IMF slope IMF slope 
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Mass- weighted age, fraction of old stellar populations (>5 Gyr) and stellar mass obtained from the derived SFH with a bimodal IMF 

shape. 



